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Abstract 

This  thesis  presents  a  comparison  of  the  longitudinal  flying 
qualities  as  predicted  by  an  analytical  computer  model,  a  ground- 
based  simulator,  and  an  airborne  simulator.  The  comparison  was 
designed  to  correlate  the  results  and  judge  whether  ground  tests 
could  forecast  airborne  results.  The  objective  was  to  determine 
and  compare  the  pilot  performance  in  the  three  cases,  and  to 
determine  and  compare  the  pilot  ratings  in  the  simulators. 
Secondary  objectives  included  an  investigation  of  the  effects  of 
different  load  factors  in  both  simulators  and  the  effects  of 
visual  and  motion  cues  in  the  NT-33A.  Identical  aircraft 
dynamics,  flight  control  characteristics,  and  tracking  tasks  were 
used  in  each  case.  The  handling  characteristics  for  short  period 
natural  frequencies  of  2,  4,  6,  8,  and  10  radians  per  second  were 
rated  using  the  Cooper-Harper  rating  scale.  In  the  analytical 
model,  pilot  performance  improved  as  the  frequency  was  increased. 
In  the  ground  simulator,  the  pilot  ratings  were  primarily  a 
function  of  how  well  he  could  track  the  preprogrammed  task. 
Performance  improved  as  frequency  was  increased.  In  the  airborne 
NT-33A  tests,  pilots  preferred  4  to  6  radians  per  second.  Lower 
frequencies  were  too  slow  and  higher  frequencies  were  too  abrupt 
and  uncomfortable  despite  better  tracking  performance.  The 
inconsistent  pilot  preference  above  6  radians  per  second  in  the 
two  simulators  is  due  to  the  absence  of  motion  cues.  The  lack  of 
correlation  at  the  higher  frequencies  indicates  that  ground  based 
simulation  cannot  entirely  replace  airborne  testing. 
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I .  INTRODUCTION 


This  thesis  presents  a  comparison  of  the  longitudinal  flying 
qualities  as  predicted  by  an  analytical  computer  model,  a  ground- 
based  simulator,  and  an  airborne  simulator.  The  comparison  was 
designed  to  correlate  the  results  and  judge  whether  ground  tests 
could  forecast  airborne  results.  Project  tests  were  aducted 
using  the  Optimal  Pilot  Single  Axis  Control  Task  (OPS  . < 
computer  program,  the  ground  based  Simulator  for  Air<  -  Flight 
Test  and  Development  ( SAFTD ) ,  and  the  variable  stabile USAF 
NT-33A  aircraft.  The  objective  of  this  evaluation  was  to 
determine  and  compare  the  pilot  performance  in  the  three  cases, 
and  to  determine  and  compare  the  pilot  ratings  in  the  simulators. 
Secondary  objectives  included  an  investigation  of  the  effects  of 
different  load  factors  in  the  two  simulators  and  the  effects  of 
visual  and  motion  cues  in  the  NT-33A 

Tight  experimental  control  was  maintained  by  programming 
identical  aircraft  dynamics  into  the  models  and  setting  optimum 
values  for  both  the  short  period  damping  and  the  control  system 
dynamics  such  as  forces,  gradients,  and  stick  displacements. 

These  values  were  obtained  from  previous  research  at  the  same 
flight  condition  which  corresponds  to  an  n/ck  of  29  g's  per 
radian.  A  stochastic  pitch  tracking  task  was  used  in  the 
computer  model  while  identical  pitch  tracking  tasks  were 
programmed  into  the  SAFTD  TV  display  and  the  NT-33A  HUD.  In  the 
OPSACT,  the  root  mean  square  of  the  tracking  error  was 


1 


determined.  In  the  two  simulators,  pilot  comments  were  recorded 
and  tracking  ability  at  short  period  natural  frequencies  of  2,  4, 
6,  8,  and  10  radians  per  second  was  rated  using  the  Cooper-Harper 
rating  scale.  Actual  tracking  error  was  also  recorded  on  strip 
charts  in  the  SAFTD  and  a  magnetic  tape  data  recorder  in  the 
NT-33A.  These  sources  provided  a  quantitative  measure  of  actual 
pilot  performance.  In  addition,  HUD  camera  film  was  also  used  to 
record  airborne  tracking  tasks . 

Computer  testing  of  the  analytical  model  took  place  on  the 
Cyber  computers  at  Wright-Patterson  Air  Force  Base,  Ohio,  and 
Edwards  Air  Force  Base,  California.  Ground  simulation  was 
conducted  in  the  SAFTD  at  Edwards  Air  Force  Base,  California, 
between  28  April  and  20  May  1983.  During  26  hours  of  SAFTD 
testing,  125  lg  test  points  and  37  2g  ^est  points  were 
accomplished.  Airborne  testing  was  conducted  at  Edwards  AFB, 
California  from  20  April  to  6  May  1983,  consisting  of  20  NT-33A 
sorties  for  a  total  of  24.4  flying  hours.  During  NT-33A  tests, 
153  lg  and  71  2g  test  points  were  accomplished  with  33  being 
flown  with  outside  visual  reference  denied  using  canopy  and  visor 
filters . 

The  Optimal  Pilot  Model  has  been  used  by  Harvey  (Ref.l)  to 
predict  the  long  term  (100  seconds)  performance  of  an  aircraft 
with  a  lead  computing  sight.  Mullen  (Ref.  5)  compared  the  short 
term  (10  seconds)  performance  prediction  of  the  0PM  to  a 
simulation  with  less  success  than  Harvey. 
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II. 


TEST  OBJECTIVE 


The  objective  of  this  evaluation  was  to  determine  and 
compare  the  pilot  pitch  tracking  performance  in  the  three  cases, 
and  to  determine  and  compare  the  pilot  ratings  in  the  simulators. 
From  this  comparison,  a  correlation  between  the  computer  model  or 
the  ground-based  simulator  and  the  inflight  simulator  would 
indicate  the  value  of  each. 

Secondary  objectives  included  determining  the  effects  of  a 
higher  load  factor  (2g)  on  pilot  opinion  in  the  two  simulators 
and  the  effects  of  visual  and  motion  cues  in  the  NT-33A.  The 
testing  at  different  load  factors  was  designed  to  test  the 
assumption  that  preferred  flying  qualities  at  lg  are  also 
preferred  at  higher  load  factors.  The  investigation  of  the 
visual  and  motion  cues  was  designed  to  determine  which  cues  are 
most  important  to  the  pilot's  ratings. 

III.  TEST  ITEM  DESCRIPTION 
THE  OPTIMAL  PILOT  MODEL 

The  optimal  pilot  model  used  in  the  analytical  portion  of 
this  study  was  the  Optimal  Pilot  Single  Axis  Control  Task 
(OPSACT).  This  computer  program  has  been  previously  used  to 
predict  pilot  opinion  of  longitudinal  tracking  tasks  (Ref.  3  and 
5). 

Based  upon  optimal  stochastic  control  theory,  the  original 
concept  for  the  model  was  developed  by  Kleinman,  Baron,  and 
Levison  (Ref.  2).  It  relies  on  the  assumption  that  a  human 
operator  with  good  training  behaves  in  an  optimal  manner  for  a 


given  control  task.  OPSACT  was  assembled  by  Enright  using  this 
concept  of  Kleinman  et  al .  (Ref.  3). 

The  model  incorporates  a  closed  loop  system  that  includes 
the  system  dynamics  of  the  aircraft  and  the  tracking  task,  the 
display  of  these  dynamics,  and  the  pilot's  perception  and  control 
actions  (Figure  1).  The  uncertainties  of  the  real  process  are 
modeled  by  stochastic  noises  which  enter  the  system  as  a  random 
target  load  factor,  errors  in  pilot  perception,  and  errors  in  the 
pilot  control  inputs. 


Figure  1.  Diagram  of  the  Optimal  Pilot  Model 


The  systems  dynamics  (Block  A)  are  a  combination  of  the  aircraft 


dynamics  and  the  pitch  tracking  dynamics.  One  of  the 
assumptions  in  the  developement  of  this  model  was  that  the  Ig 
flying  qualities  of  an  aircraft  are  not  significantly  altered 
when  the  load  factor  is  increased.  Appendix  A  lists  the  results 
of  an  investigation  that  showed  this  to  be  true  as  the  frequency 
response  of  the  NT-33A  aircraft  was  not  change  significantly  as 
the  load  factor  was  increased. 

The  tracking  task  used  in  this  model  was  a  stocastic  pitch 
tracking  task  designed  to  simulate  the  pitch  tracking  task  that 
the  pilots  flew  in  the  simulators.  It  is  a  zero  mean  pitch 
command  with  a  deviation  and  time  constant  similar  to  the  task 
programmed  into  the  SAFTD  display  and  the  NT-33A  HUD. 

The  state  equations  of  motion  were  formed  using  the 
following  state  variables: 

xx  -  $ 

X2  =  q 

x3  =  <* 

and  u  =  £  =  -K£F 

uc  f  s 

The  equations  of  motion  and  target  dynamics  are  combined  into 
matrix  notation 

x  =  Ax ( t )  +  Bu(t)  +  W(t) 


x4  =  © 

=  ©t 

x6  =  dummy  noise  state , 
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in  which  the  "A"  matrix  is 


M  +M42t 
Uo 


M^+Mi.Z. 


-KfFs 


For  300  KIAS  and  W  =  2  radians  per  second,  the  A  matrix  is 

sp 


-22-28  -0.994  -2-204 


-0.105 


-1.81 


Also,  B  =  col  [0000001] 


-0.40 


0.114 


-1.339 


and  W(t)  =  col  [0  0  0  0  0  w  0] 

where  w  =  white  Gaussian  driving  noise. 

The  pilot  sees  the  pitch  tracking  error  (Block  B),  £.  which 
is  e-Sc,  and  the  rate  of  change  of  the  error,  £.  The  model  of 


this  is 


y(t)  =  H  x ( t ) 
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with  y(t)  =  col  [6  €  ] 

and  the  output  or  observation  matrix  H  is 

00  0  1  -1  o  o" 

_  0  1  0  0  l/t  -k  0 

The  time  delays  that  are  associated  with  the  relaying  and 
The  time  delays  that  are  associated  with  the  relaying  and 
processing  of  visual  images  in  the  brain  are  modeled  by  lumping 
them  into  a  single  equivalent  time  delay  (Block  C) .  This  delay, 
'T ,  is  nominally  0.2  jt0.05  seconds. 

A  Kalman  estimator  (Block  D)  works  in  series  with  a  least 
mean  square  predictor  (Block  E)  to  yield  a  "best  estimate"  of  the 
system  state  vector  from  a  time  delayed,  noisy  observation 
vector.  The  estimate  is  then  weighted  by  a  set  of  optimal 
feedback  gains  which  are  determined  from  the  solution  of  an 
optimal  control  problem.  This  weighted  estimate  is  used  to 
produce  a  scalar  command  control  input  (Block  F),  uc(t),  which  is 
a  longitudinal  stick  force. 

The  neuromotor  dynamics  are  modeled  by  including  a  first 
order  lag  (Block  G) .  This  lag  accounts  for  the  pilot’s  inability 
or  reticence  to  make  rapid  or  excessive  control  inputs. 

Controller  remnant  accounts  for  the  inherent  random  errors 
associated  with  the  preception  of  displayed  variables  and  in  the 
pilot  control  inputs. 

Since  the  model  is  linear,  these  errors  are  consolidated  as 

observation  noises,  v  ,  and  motor  noise,  v  .  These  v^  and  v  are 

V  u  y  u 

assumed  independent,  zero  mean  Gaussian  noises  with  enough 
bandwidth  to  be  considered  white  noise  processes. 
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A  detailed  description  of  OPSACT  is  contained  in  References 
3  and  5 . 

THE  GROUND  SIMULATOR  ( SAFTD) 

The  Simulator  for  Aircraft  Flight  Test  and  Development 
(SAFTD)  is  a  real-time  simulator  laboratory  used  to  support 
flight  test  programs  at  the  Air  Force  Flight  Test  Center.  The 
Test  Pilot  School's  simulator  model  was  developed  from  an 
existing  F-16  model.  Two  digital  computers  arranged  in  a 
distributive  processing  network  were  used  with  a  generic  fixed 
base  fighter  cockpit  with  a  moveable  center  stick  as  shown  in 
Figure  2.  The  stick  was  connected  to  a  McFadden  artificial  feel 
system  which  generated  the  stick  dynamics.  Simulator  operation 
was  controlled  using  an  interactive  alphanumeric  cathode  ray  tube 
(CRT)  display  and  data  entry  keyboard.  Two  synchronous 
processors  used  aerodynamic  coefficient  lookup  tables  while 
integrating  the  appropriate  equations  of  motion.  Thus  the 
simulator  modeled  the  aircraft  for  all  values  of  airspeed  and 
angle  of  attack  (AO A) .  A  series  of  test  pulses  identical  to 
those  used  in  the  NT-33A  were  fed  to  a  pitch  steering  bar 
displayed  on  a  25  inch  raster  scan  TV.  The  25  inch  TV  was 
partially  covered  to  match  the  dimensions  of  the  NT-33A  HUD.  The 
maximum  time  delay  between  the  stick  input  and  update  of  the  TV 
display  was  three  frames  or  48  milliseconds.  This  is  well  below 
the  pilot  perceptual  threshold  of  approximately  100  milliseconds 
(Ref.  18).  A  time  history  of  eight  key  flight  parameters  were 
displayed  on  strip  charts  (Appendix  D). 
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Figure  2.  SAFTD  Layout 


THE  NT-33A 

The  test  aircraft,  NT-33A  Serial  Number  51-4120,  is  a  T-33A 
modified  with  a  variable  stability  system  by  CALSPAN  Corporation, 
Buffalo,  New  York.  This  variable  stability  system  uses  response 
feedback  to  modify  the  static  and  dynamic  responses  of  the  basic 
NT-33A  by  commanding  control  surface  positions  through  full 
authority  electrohydraulic  servos.  System  components  include  a 
programmable  analog  computer,  associated  aircraft  response 
sensors,  control  surface  servos,  and  an  electrohydraulic 
force-feel  system.  The  rear  seat  pilot  can  vary  the  computer 
gains  through  controls  located  in  the  rear  cockpit  allowing 


changes  in  airplane  dynamics  and  control  system  characteristics 
in  flight.  Stick  friction  and  breakout  forces,  stick  deflection 
gradients,  and  stick  force  gradients  were  achieved  through  an 
electrohydraulic  system.  A  schematic  of  the  aircraft  flight 
control  system  is  presented  in  Figure  3. 


Figure  3.  NT-33A  Control  System  Layout 

An  on-board  Leach  MTR  3200  magnetic  tape  recording  system 
was  utilized  to  record  aircraft  flight  conditions,  flight  control 
positions,  and  the  tracking  task  parameters.  A  front  cockpit 
AVQ-7  HUD  displayed  the  tracking  task.  The  system  incorporated  a 
general  purpose  digital  computer  and  programmable  display 
generator.  This  system  has  a  processing  time  delay  of  32 
milliseconds.  The  pitch  tracking  task  that  appeared  on  the  HUD 
was  preprogrammed  on  magnetic  tape.  A  16mm  HUD  camera  was  used 


\ 
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IV.  TEST  METHODS  AND  CONDITIONS 

OVERALL 

An  important  element  in  the  achievement  of  valid  results  was 
the  matching  of  the  dynamics  of  the  various  tests.  The 
analytical  model  used  the  short  period  approximation  to  model  the 
aircraft.  Baseline  testing  was  accomplished  in  the  SAFTD  and 
NT-33A  to  verify  that  the  configuration  of  each  matched  the 
dynamics  of  the  short  period  approximation.  This  was 
accomplished  by  analog  matching  open  loop  strip  chart  data  with 
computer  generated  response  plots.  Figure  5  shows  an  example  of 
the  analog  matching. 


Computer  Trace 

Response  to  Pitch 


SAFTD 
Doublet,  ^Sp 


NT-33A 
6  Rac/Sec,  ^  =  0.7 


Figure  5.  Sample  Analog  Match 


Also,  the  breakout  force,  friction  force,  stick  deflection 
gradient,  stick  force  gradient,  and  the  stick  damping  in  the 
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simulators  were  compared  and  verified  using  force  gauges  and 
strip  charts. 

ANALYTICAL  MODEL 

Tests  with  the  analytical  model  were  conducted  using 
aerodynamic  derivatives  from  the  NT-33A  at  the  flight  conditions 
as  the  inflight  tests.  These  conditions  corresponded  to  300  KIAS 
at  13,000  feet  MSL .  This  gives  a  normal  load  factor  per  angle  of 
attack  (n/cx)  of  29  g’s  per  radian. 

The  short  period  natural  frequency  was  varied  by  changing  M^ 
with  being  changed  to  maintain  a  constant  0.7  damping  ratio. 

The  elevator  linkage  gain  was  varied  to  maintain  a  constant  stick 
force  per  g.  This  corresponded  to  the  technique  used  in  the 
SAFTD  and  the  NT-33A. 

SIMULATION 

Data  collection  was  accomplished  by  performing  identical  HUD 
pitch  tracking  tasks  at  the  same  flight  condition  in  both  the 
SAFTD  and  NT-33A  simulators  at  specific  short  period  natural 
frequencies.  During  each  90-second  task,  only  one  frequency  was 
tested  and  tracking  ability  was  then  rated  using  the 
Cooper-Harper  rating  scale  (Figure  6)  and  pilot  comments. 


ADEQUACY  FOR  SELECTED  TASK  OR 
REQUIRED  OPERATION 


AIRCRAFT 

CHARACTERISTICS 


DEMANDS  ON  THE  PILOT 
IN  SELECTED  TASK  OR  REQUIRED  OPERATION 


PILOT 

RATING 


Excellent 

Highly  desirable 

Pilot  compensation  not  a  factor  lor 
desired  pertor monte 

n 

Good 

Negligible  deficiencies 

Pilot  compensation  not  a  factor  lor 
desired  performance 

Fa«r  —  Some  mildly 
unpleasant  deficiencies 

Minimal  pilot  compensation  required  lor 
desired  pertor monce 

Minor  but  annoying 
deiiciencics 

Desired  performance  requires  moderate 
pilot  compensation 

pTl 

Moderately  obieciionable 
dehciencies 

Adequate  performance  requires 
considerable  pilot  compensation 

Very  objectionable  bul 
tolerable  deficiencies 

Adequate  performance  requires  extensive 
pilot  compensation 

Major  deficiencies 

Adequate  performance  not  attainable  with 
maximum  tolerable  pilot  compensation 
ControHabil'ty  not  m  question 

n 

Major  deficiencies 

Considerable  pilot  compensation  is  required 
for  consol 

Maior  deficiencies 

Intense  pilot  compensation  is  required  to 
retain  control 

Ijjjl 

Maior  deficiencies 

Control  will  be  tost  during  some  portion  of 
required  operation 

Figure  6.  The  Cooper-Harper  Rating  Scale 

The  pitch  tracking  tasks  consisted  of  a  90-second  series  of 
steps  and  ramps  which  commanded  up  to  +  4  degrees  from  the 
initial  level  flight  reference.  The  two  different  pulse  trains, 
presented  in  Figure  7,  were  used  to  prevent  pilots  from 
memorizing  the  task.  Both  tasks  were  designed  by  CALSPAN  to 
excite  longitudinal  dynamic  modes  that  were  important  for 
successful  air-to-air  pitch  tracking. 
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Figure  7.  Preprogrammed  Pitch  Tracking  Tasks 


The  pitch  command  bar  was  programmed  to  move  in  pitch 
according  to  the  tracking  sequence.  The  SAFTD  task  was  displayed 
an  a  television  screen  while  the  NT-33A  display  was  on  a  HUD 
combining  glass.  The  project  pilot  attempted  to  track  the  moving 
command  bar  with  the  top  of  the  "W"  (waterline  marker)  in  both 
simulators.  The  "W"  was  the  same  size  in  each  simulator  (5.3 
milliradians) .  A  moveable  centerstick  controller  was  used  for 
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all  tests  in  both  simulators.  Desired  performance  in  both 
simulators  was  defined  as  acquisition  of  the  moving  pitch  command 
bar  with  the  top  of  the  "W"  within  one  second  and  then  maintain 
pitch  attitude  within  a  "W"  width  for  the  duration  of  the  test 
pulse.  Adequate  performance  was  acquisition  within  two  seconds 
and  then  maintain  pitch  attitude  within  two  "W"  widths  of  the 
pitch  command  bar  for  the  duration  of  the  pulse. 

The  tracking  task  was  repeated  in  each  simulator  varying 
only  the  short  period  natural  frequency  (WSp) ,  which  was  varied 
in  even  increments  from  2  to  12  radians  per  second  in  the  SAFTD 
and  from  2  to  10  radians  per  second  in  the  NT-33A.  Twelve 
radians  per  second  could  not  be  tested  in  the  NT-33A  because  of 
structural  mode  excitation.  At  this  frequency,  an  unstable  AOA 
feedback  to  the  elevator  caused  the  variable  stability  system  to 
exceed  its  design  limits  and  automatically  disengage.  All  other 
aircraft  and  control  system  variables  were  selected  to  provide 
good  flying  qualities.  These  optimized  values  were  obtained  from 
previous  testing  at  an  n/c*  of 

29  g's  per  radian  and  are  listed  in  Table  1  (Ref  12).  The 
lateral  directional  characteristics  were  chosen  to  provide  good 
control  harmony. 

The  ratio  of  n/crt  was  maintained  at  29  g's  per  radian  in  both 
simulators  by  controlling  altitude  and  airspeed  at  13,000  +  1,000 
feet  pressure  altitude  and  300  +  15  KIAS  respectively.  The  SAFTD 
was  programmed  to  simulate  the  NT-33A  at  this  flight  condition. 
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Table  1 

SELECTED  DYNAMIC  PARAMETERS 
(13,000  ft  PA,  300  KIAS) 


PARAMETER 

VALUE 

Shoct  period  natural  frequency 

2-12 

(rad/sec) 

Short  period  damping 

0.70 

Normal  g  per  alpha 

29.00 

(g/rad) 

Dutch 'roll  natural  frequency 

3.20 

(rad/sec) 

Dutch  roll  damping 

0.35 

Phi  to  beta  ratio 

2.00 

Roll  mode  time  constant 

0.35 

(sec) 

Phugoid  natural  frequency 

0.09 

(rad/sec) 

Phugoid  damping 

0.05 

Longitudinal  stick  force  per  g 

6.50 

( lb/g ) 

Longitudinal  stick  force  per  inch 

8.00 

(lb/in) 

Test  points  were  selected  to  allow  equal  pilot  exposure  to 

each  configuration.  The  value  of  OJ  was  unknown  to  the  project 

sp 

pilot  and  the  sequencing  order  of  ^Sp ' s  tested  was  determined 
with  a  random  number  generator. 

Limited  testing  was  done  in  both  simulators  at  a  nominal 
2g's  by  trimming  for  lg  level  flight  and  then  maintaining  a 
nominal  60  degree  banked  turn  throughout  the  tracking  sequence. 

The  effects  of  outside  visual  cues  on  pilot  ratings  in  the 
NT-33A  were  investigated  by  denying  outside  references. 

Temporary  restricted  vision  was  accomplished  with  a  blue  helmet 
visor  and  an  amber  NT-33A  cockpit  screen.  This  combination  made 
the  HUD  clearly  visible  while  making  the  scene  beyond  the  HUD 
appear  as  darkness. 


V.  TEST  RESULTS 


Test  results  are  presented  in  three  parts:  the  OPSACT 
results,  the  SAFTD  results,  the  NT-33A  results,  and  a  comparison 
of  the  three. 

In  the  OPSACT  model,  the  performance  of  the  "pilot"  got 
progressively  better  as  the  frequency  was  increased  (Figure  8). 
The  Qy  matrix  is  the  relative  weight  of  the  tracking  error  and 
the  rate  of  change  of  the  tracking  error. 


poor  and  gave  successively  better  ratings  as  short  period 


frequency  was  increased.  Representative  pilot  comments  for 
tested  short  period  frequencies  are  summarized  in  Table  2. 


Table  2 

SAFTD  GROUND  SIMULATOR  PILOT  COMMENTS 

U)  sp 

COMMENTS 

2 

Very  slow,  poor  predictability,  large 
overshoots,  sluggish 

4 

Slightly  slow,  fairly  predictable, 
moderate  overshoots,  large  inputs 

6 

Quick,  moderate  predictability,  small 
overshoots,  lots  of  compensation 

8 

Quick,  good  predictability,  small 
oscillations 

10 

Really  quick,  good  predictability, 
small  overshoot,  hardly  any  PIO 

12 

Very  quick,  good  predictability,  small 
overshoot,  pretty  comfortable,  no  PIO 
problem 

In  addition  to  the  pilots'  word  picture,  graphs  of 
Cooper-Harper  ratings  versus  short  period  frequency  are  plotted 
in  Figures  9A  and  9B. 
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Figure  9B.  Cooper-Harper  Ranges  for  SAFTD  at  lg 


Figure  9A  depicts  statistical  confidence  intervals  showing 
the  predominant  Cooper-Harper  ratings  for  each  tested  frequency. 
It  is  not  possible  to  plot  exact  statistical  intervals  such  as  90 
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or  95  percent  for  each  frequency  as  the  Cooper-Harper  ratings  are 
discrete  values.  The  confidence  interval  was  determined  by  the 
formula : 

S-i  p 

Pr(Yr  <*«•  <Ys  )  =  (  )  ( 1/2  )n 

This  states  that  there  is  a  fixed  probability  that  OL  the  median 
Cooper  Harper  rating!  will  be  between  a  minimum  and  maximum 
rating,  Y  and  Y  .  If  the  ratings  are  rank  ordered,  r  is  the 
position  of  the  first  occurrence  of  the  lower  rating  and  s  is  the 
last  occurrence  of  the  higher  rating.  The  relationship  applies 
regardless  of  the  distribution  of  the  sample  or  the  value  of  the 
median.  Figure  9B  shows  the  absolute  range  of  Cooper-Harper 
ratings  for  each  tested  frequency  as  well  as  the  median  and  mean 
value.  Figures  10A  and  10B  present  the  2g  test  results  for  the 
SAFTD  in  the  same  format . 


OMEGA  <W>.  NATURAL  FREQUENCY  (rad/sec) 

Figure  10A.  Confidence  Intervals  for  SAFTD  at  2g 


21 


Figure  10B.  Cooper-Harper  Ranges  for  SAFTD  at  2g 


The  SAFTD  results  show  an  increasing  pilot  preference  for 
higher  frequencies.  SAFTD  rating  confidence  intervals  show  a 
negligible  difference  between  lg  and  2g  testing  while  the  mean 
ratings  are  worse  at  2g .  The  preferred  frequencies  are  the  same 
at  lg  and  2g,  supporting  the  assumption  that  flying  qualities 
that  get  the  better  ratings  at  low  load  factors  will  get  the 
better  ratings  at  higher  load  factors. 

Figure  11  shows  how  the  SAFTD  pilot  performance  inproved  as 
the  short  period  natural  frequency  was  increased. 
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Figure  11.  SAFTD  Tracking  Errror  Versus  Frequency 


The  NT-33A  in-flight  results  follow  other  flight  test 
results.  The  pilots  disliked  the  very  low  frequencies,  liked  the 
predictable  medium  response  area,  and  disliked  the  uncomfortably 
abrupt  high  frequencies.  As  in  the  SAFTD,  there  was  no 
significant  difference  in  pilot  opinion  and  Cooper-Harper  ratings 
between  lg  and  2g  tests,  with  the  numerical  differences  primarily 
attributed  to  sample  size.  Here,  there  was  not  the  difference  in 
mean  ratings  that  occurred  in  the  SAFTD.  Testing  with  no  outside 
visual  references  had  a  negligible  effect  on  ratings.  The  pilots 
focused  on  the  HUD  during  the  tracking  task  in  either  situation. 


eliminating  background  visual  cues  as  a  major  factor  in  pilot 
ratings . 

The  NT-33A  lg  results  are  presented  in  Table  3  and  Figures 
12A  and  12B.  Table  3  summarizes  the  pilot  opinions  concerning 
response,  predictability,  PIO  tendency,  and  work  load.  Figure 
12A  presents  the  statistical  confidence  interval  on  the  median 
Cooper-Harper  rating  and  Figure  12B  shows  the  absolute  range  of 
ratings  given  to  each  tested  frequency  as  well  as  the  sample 
median  and  mathematical  mean. 


Table  3 

NT-33A  PILOT  COMMENTS 

U>  Sp 

COMMENTS 

2 

Sluggish,  poor  predictabil ity , 
large  overshoots,  a  lot  of 
compensation 

4 

Good  response ,  good  predictability, 
comfortable,  minimal  compensation 

6 

Good  response,  not  abrupt,  good 
predictability,  feels  good 

8 

Quick  response,  abrupt,  less 
predictable,  a  little  uncomfortable, 
some  PIO  tendency 

18 

Jerky,  too  jumpy,  increased  work 
load,  lacge  compensation, 
uncomfortable,  bobbles,  PIO  prone 

24 
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Figure  12A.  Confidence  Intervals  for  NT-33A  at  lg 


Figure  12B.  Cooper-Harper  Ranges  for  NT-33A  at  lg 


The  2g  test  results  for  the  NT-33A  are  presented  in  Figures 
13A  and  13B.  Like  the  SAFTD  results,  there  was  a  negligible 
difference  between  lg  and  2g  ratings. 
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Figure  13A.  Confidence  Intervals  for  NT-33A  at  2g 


2  4  6  8  •  10  12 

OMEGA  <M^,  NATURAL  FREQUENCY  (rad/sec) 


Figure  1 3B •  Cooper-Harper  Ranges  for  NT-33A  at  2g 

Figure  14  is  a  plot  of  pilot  performance  versus  short  period 
natural  frequency  in  the  NT-33A.  rms  error  was  not  available 
because  of  format  incompatabilities  between  the  NT-33A  data 
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acquisition  system  and  the  Test  Pilot  School  computer.  Error 
counts  were  manually  produced  from  strip  charts  by  determining 
the  area  under  the  curve  which  was  the  error  between  the  pitch 
command  and  the  actual  pitch  angle.  This  provides  the  same 
relative  indication  of  performance  versus  frequency  as  RMS  error. 


103  0  z  4  6  8  10  12 

NAT  FREQ  (RRD/SEC  >  , 

Figure  14.  NT-33A  Error  Count  Versus  Frequency 

The  project  objective  was  to  compare  the  OPSACT  results,  the 
SAFTD  results , and  the  NT-33  in-flight  results  and  see  if  ground 
tests  could  forecast  airborne  opinions.  Tight  experimental 
control  was  maintained  so  that  any  correlation  between  the  ground 
and  in-flight  simulators  was  not  biased  or  circumstantial. 
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VI.  ANALYSIS 


This  section  offers  an  explanation  for  the  considerable 
difference  at  the  higher  frequencies  between  the  ground-based 
studies  and  the  airborne  simulation  despite  each  case  having 
identical  equations  of  motion  and  control  dynamics. 

Performance  error  or  root  mean  square  (RMS)  error  was 
measured  as  the  difference  between  the  preprogrammed  pitch  angle 
commanded  and  the  actual  pitch  angle  which  the  pilot  controlled. 

A  small  RMS  error  figure  signifies  that  the  command  bar  was 
quickly  and  accurately  acquired  while  a  large  number  implies  slow 
response  and/or  tracking  difficulty.  Charting  performance  was 
not  a  measure  of  the  "best  pilot",  but  rather  provided  insight  as 
to  why  a  pilot  preferred  certain  frequencies  and  could  indicate 
whether  a  pilot  performed  consistently  in  different  runs  testing 
the  same  frequency. 

The  OPSACT  allows  weighting  of  the  pilot's  control  inputs 
due  to  tracking  error  and  rate  of  change  of  tracking  error.  This 
study  found  that  any  weight  greater  than  one  percent  on  the  rate 
cause  the  pilot  to  prefer  the  lower  frequency  of  two  radians  per 
second.  If  the  weighting  was  such  that  that  the  OPSACT 
performance  at  two  radians  per  second  paralleled  the  SAFTD  and 
NT-33A,  the  OPSACT  results  followed  the  SAFTD. 

Since  the  SAFTD  could  not  p»-'. '  ide  "feel"  cues,  the 
performance  criterion  dominated  p^ot  opinion  at  all  frequencies 
during  the  ground  test  phase.  In  the  SAFTD,  Cooper-Harper 
ratings  were  directly  proportional  to  RMS  performance  which 
improved  with  successively  higher  frequency.  Figure  16  depicts 
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the  pilots'  performance  versus  the  Cooper-Harper  ratings 
assigned.  This  evidences  linear  behavior  with  an  80  percent 
correlation.  Figure  11  depicts  the  RMS  error  versus  short  period 
frequency  as  a  second  order  behavior.  The  graph  indicates  that 
performance  improves  as  frequency  increases.  A  second  order  fit 
had  a  74  percent  statistical  correlation.  This  suggests  the 
hypothesis  that  performance  is  the  primary  variable  affecting 
pilot  ratings  in  the  ground  simulator.  Unlike  the  NT-33A,  there 
are  no  uncomfortable  motion  cues  at  the  higher  frequencies. 
Therefore,  Cooper-Harper  ratings  in  the  SAFTD  tend  to  improve  as 
frequency  increases. 
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Figure  16.  SAFTD  Tracking  Error  Versus  Cooper-Harper 
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In  the  NT-33A  tests,  the  pilot  opinions  were  driven  by  three 
factors:  tracking  accuracy  (performance),  comfort,  and  pilot 
compensation.  These  factors  changed  in  relative  importance  as 
the  frequency  changed.  The  pilots'  comments  indicated  that 
performance  or  ability  to  accurately  track  the  pitch  command  bar 
was  a  primary  variable  in  determining  the  flying  quality  ratings 
for  low  and  medium  frequencies.  At  higher  frequencies,  however, 
"jerkiness"  or  uncomfortably  abrupt  motions  forced  the  pilot  to 
make  large  compensations  and  drove  the  ratings  to  worse  levels. 
Plots  of  performance  versus  Cooper-Harper  ratings  and  short 
period  frequency  in  the  NT-33A  appear  in  Figures  17  and  14. 
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Figure  17.  NT-33A  Error  Count  Versus  Cooper-Harper 


Unlike  the  SAFTD  results,  where  a  direct  linear  relationship 
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between  Cooper-Harper  ratings  and  performance  was  established, 
NT-33A  results  did  not  correlate  well  as  other  factors  affected 
the  ratings.  One  of  these  factors  was  g  onset  rate,  which  was 
rapid  enough  at  high  frequencies  to  cause  problems  in  maintaining 
steady  stick  inputs,  and  in  some  cases,  made  focusing  on  the  HUD 
difficult.  Despite  the  abrupt  motion  cues,  one  pilot  showed  the 
same  increase  in  performance  with  higher  frequency  that  the 
ground  simulator  results  suggested.  Figure  18  presents  this 
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Figure  18. 


NATURAL  FREQUENCY  <RR0/5EC> 

Pilot  #2  Error  Count  Versus  Frequency 


Hand  reduced  tracking  error  data  indicated  two  different 
behavior  patterns  at  10  radians  per  second,  the  highest  tested 
frequency  in  the  NT-33A.  Either  the  pilot  aggressively  tracked 
the  commanded  task  and  ignored  the  very  abrupt  motions  or  the 
aggressive  tracking  was  tempered  to  maintain  a  comfort  level. 
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The  resulting  error  counts  reflect  this  different  level  of 
aggressiveness.  Therefore,  it  is  possible  that  pilot  performance 
could  be  improved  with  frequency  if  the  effects  of  the  other 
factors  could  be  mitigated.  This  "other  factor"  hypothesis  is 
amply  evidenced  by  the  low  correlation  between  RMS  versus 
Cooper-Harper  plot  in  Figure  17. 

The  Cooper-Harper  rating  system  defines  major  thresholds 
based  on  a  combination  of  performance,  work  load,  and  pilot 
compensation.  Thus,  performance  is  not  an  overriding  factor  in 
pilot  opinion.  Further  research  on  analytical  modeling  and 
ground  simulation  needs  to  incorporate  items  such  as  pilot 
workload  and  pilot  compensation  to  keep  the  g  onset  rate 
comfortable  into  the  prediction. 

There  was  no  significant  difference  between  lg  and  2g 
results  during  project  tests.  These  results  are  consistent  with 
the  analysis  shown  in  Appendix  A.  Therefore,  the  Military 
Specification  contention  that  lg  predictions  characterize  higher 
load  factors  cannot  be  refuted.  Since  the  project  was  primarily 
aimed  at  a  lg  comparison  of  simulators,  2g  was  only  tested  for  a 
statistically  significant  sample  size.  The  loaded  testing  should 
be  continued  at  higher  load  factors  and  in  different  aircraft. 
Testing  in  the  NT-33A  at  the  300  KIAS,  13,000  feet  MSL  flight 
condition  at  no  more  than  2g  allowed  he  aircraft  to  remain  in 
the  linear  portion  of  the  lift  curve.  Testing  at  higher  g  in  the 
non-linear  portion  of  the  curve  or  testing  aircraft  with 
different  lift  slopes  may  produce  a  different  conclusion  than 
this  limited  evaluation. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  computer  pilot  model  and  ground  based  simulator  modeled 
the  airborne  pilot  performance  trends  well.  However,  performance 
is  not  the  only  criteria  that  a  pilot  uses  to  rate  the  flying 
qualities  of  an  aircraft.  The  pilot  considers  other  factors  such 
as  g  onset  rate  when  he  forms  his  opinion.  Therefore,  any  model 
which  uses  performance  as  a  dominate  criteria  will  not  generally 
be  successful.  The  model  can,  however,  be  used  in  parallel  with 
the  airborne  simulator  or  test  article  to  reduce  airborne  test. 

Further  work  with  the  weighting  matrix  in  the  analytical 
model  might  produce  a  combination  that  would  more  closely 
parallel  the  airborne  tests.  One  method  of  finding  this 
weighting  would  be  to  generate  a  time  history  of  the  pilot  model. 
By  comparing  the  time  history  to  a  time  history  of  the 
simulators,  the  behavior  of  the  model  could  be  adjusted  to  be 
similar  to  actual  pilots. 

The  assumption  that  the  preferred  lg  flying  qualities  will 
also  be  preferred  at  2g  and  higher  cannot  be  refuted  by  this 
study,  and,  at  least  for  longitudinal  dynamics,  is  supported  by 
the  analysis  in  Appendix  A.  More  testing  at  higher  load  factors 
is  recommended . 

For  this  type  of  task,  which  has  the  high  fidelity  of  a  HUD, 
denying  outside  visual  references  did  not  have  a  significant 
effect  on  pilot  performance  or  ratings.  This  means  that  the 
difference  between  the  SAFTD  and  NT-33A  were  dominated  by  the 


\ 
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effects  of  motion.  Therefore,  using  a  HUD  type  task  with  high 
fidelity  reduces  the  number  of  variables  between  the  studies. 
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APPENDIX  A 

ANALYSIS  OF  THE  EFFECTS  OF  LOAD  FACTOR 
ON  FREQUENCY  RESPONSE  OF  THE  NT-33A 


INTRODUCTION 


When  doing  aircraft  design  studies,  the  decoupling  of  the 
longitudinal  and  lateral  axes  makes  the  problem  tractable. 
However,  many  practical  flight  conditions  require  banked,  turning 
flight  which  causes  coupling  between  the  axes  of  motion.  By 
investigating  frequency  response  of  a  NT-33A  in  several 
steady,  level  turns,  this  appendix  determines  the  validity  of 
using  decoupled  equations  to  study  flight  conditions  where  the 
axes  are  coupled. 

Eight  interdependent  equations  of  motion  are  required  to 
model  three  dimensional  flight.  This  model  is  based  on  three 
force  equations,  three  moment  equations,  and  two  kinematic 
relationships  which  were  all  linearized  about  an  equilibrium 
condition . 

The  data  for  this  analysis  were  extracted  from  NASA  CR-2144 
(Ref.  6).  The  NASA  data  are  listed  in  two  formats,  graphs  and 
tabular  charts.  The  graphical  data  were  non-dimensional  and  in 
stability  axes.  The  tabular  data  were  dimensional  and  in  body 
axes.  To  accomplish  the  analysis,  all  data  were  transformed  into 
dimensional  form  and  set  in  the  stability  axis  system  appropriate 
for  straight  and  level  (lg)  flight.  The  change  in  angle  of 
attack  required  to  maintain  each  load  factor  was  then  determined. 
Finally,  the  data  were  transformed  to  the  sets  of  stability  axes 
corresponding  to  each  load  factor  and  the  frequency  response  for 
the  appropriate  control  input  was  determined  using  TOTAL,  a 
computer  program  for  analyzing  linear  dynamical  systems  (Ref.  7). 


VARIABLES  OF  INTEREST 


The  important  variables  in  a  handling  qualities  study  are 

those  that  are  easily  preceived  and  controlled  by  the  pilot. 

Rate  variables  such  as  q,  the  pitching  rate,  and  p,  the  rolling 

rate,  are  not  affected  by  the  orientation  of  the  coordinate 

system.  Thus,  q  and  p  were  chosen  as  the  variables  of  interest. 

Theta  was  chosen  so  it  could  be  compared  to  q. 

-3  2 

The  frequency  range  of  10  to  10  radians  per  second  was 
chosen  as  this  encompassed  the  response  characteristics  of  a 
typical  pilot. 

EQUATIONS  OF  MOTION 

Assuming  a  flat,  nonrotating  earth,  a  rigid,  symmetric 
vehicle,  and  a  stationary  atmosphere  makes  the  following 
equations  applicable  to  the  problem. 

Force  Equations:  X  -  mg  sin  0  =  m  ( u  +  qw  -  rv) 

Y  +  mg  cos  8  sin  <J>  =  m  (v  +  ru  -  pw) 

Z  +  mg  cos  9  cos  $  =  m  (w  +  pv  -qu) 


Moment  Equations:  L  = 

V  - 

Izx 

(t  + 

pq) 

-  ( 

Iy 

V 

qr 

M  = 

V*  - 

Jzx 

u2 

-  P2> 

- 

(Iz 

-  Ix) 

rp 

N  = 

V  - 

*zx 

(p  - 

qr) 

-  ( 

rx  - 

V 

pq 

Kinematic  Relations: 

‘<j>  *  p  + 

q 

sin  (j, 

tan 

8  + 

r  cos  (j, 

tan 

B  =  q  cos 

♦  ~  r 

sin 

<t> 

*  =  (q 

cos 

+ 

r  cos 

♦> 

sec 

e 

Note 

:  The 

’'J' 

equa  1 

:ion 

is 

independ 

lent 

other  equations  so  it  is  not  required  for  this 
study. 
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These  equations  were  linearized  about  a  trim  condition 
representing  a  steady,  level  turn.  The  following  substitutions 
were  made: 

1.  ©e  =  0  by  definition  of  a  steady,  level  turn. 

2.  cos  /3Q  =  1  and  sin  /?e  =  0  as  the  largest  /?e  for  the  load 
factors  investigated  was  /3e  =  0.00253  radians. 

3.  p  =0  when  ©  =  0. 

rG  “ 


The  change  in  angle  of  attack  for  each  flight  condition  was 
determined  by  (Ref.  8:426): 


The  change  in  angle  of  sideslip  was  checked  using: 


The  equilibrium  values  of  p,  q,  and  r  ace  (Ref  3:424) 


The  computer  programs  which  solve  these  sets  of  equations 
are  listed  in  Appendix  B. 


1 


The  resulting  equations  are: 


[m]  AV  =  [Tv  cos  ote  -  Dv]  AV  +  [mReVe]  AB 

+  t"Te  sin  ote  +  Le  -  Da  -  mqeVe]  Aa 

+  [-mg]  A0  +  [D  .  ]  A6  e 
Se 

[mVe ]  aF  =  [mrel  AV  +  [Y  1  AB  +  [Ypl  AP 

+  l-mVe  +  Yr]  ac  +  [mg  cos  $  g]  a<J>  +  [Y  J  A6 


[mVe  +  L  .]  Aa  =  t>mqe  "  Tv  sin~  te 


"  L„]  AV 


+  [-Te  cos  a(.e 


"  De  "  La]  Act  +  tmVe  _  Lq]  Aq 


+  [-mg  sin  <j>e]  a$  +  [L  A6  e 
Ux]  A?  +  l-I2X]  A?  =  [Lgl  AS  +  t^p  +  Izxqe^  AP 

+  Kiy  -  Iz)  tei  Aq 
+  CLr  +  (iyy  -  izz)  qe3  Ar  +  [l*]  a£ 

[Iv]  A^  +  [-Ml  Aa"  =  [M„]  AV  +  [M  ]  Aa 
a  a  v  a 

+  mz  -  Ix)  reJ  AP  +  [Mq]  Aq 
+  [2IZxre]  AC  +  tM5e]  A6e 

Hz]  A?"  +  t"Izx]  AP  =  IN  ]  AB  +  [Np  +  (Ix  -  Iy)  q0]  AP 

+  l-IzxCe]  Aq  +  [Nr  -  Izxqe]  At 
+  [Nfi]  AS 

a|"  =  AP  +  tqe  sin  <f>  e  +  re  cos<f>  e1  Ae 


A0 


[cos  a  ]  ^  +  t-sin  <|>  ]  Ar  +  [-q  sin  <j>  -  r  cos  <J>  ]  A<J> 
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DATA  MAN  I PULAT ION 

The  following  algorithm  was  computerized  for  the  transfer  of 
data  from  its  raw  form  to  that  required  by  TOTAL: 


FLOW  CHART  FOR  DATA  PROCESSING 


/  Enter  constants 

and  Non-dimensional  Quantities/ 


Convert  Inertia  Terms  to  Stability  Axes 
Dimensionalize  All  Non-dimensional  Data 


Enter  Dimensional  Data, 


FREQUENCY  RESPONSE 


* 


The  program  TOTAL  takes  equations  of  the  form 
x  =  [A]  x  +  [B]  u 
y  -  CC]  x 

and  generates  the  transfer  functions  between  user  designated 
inputs  and  outputs.  From  these  transfer  functions,  the  program 
is  able  to  generate  Bode  plots. 

RESULTS 

The  Bode  plots  in  Appendix  B  show  the  frequency  response  of 
the  NT-33A. 

In  each  case,  the  ©/£,,  response  is  of  the  same  general  form 
with  a  spike  at  the  damped  frequency  of  the  phugoid.  This  occurs 
at  higher  frequency  as  the  loa^  factor  is  increased  (Figure 
B1-B5).  The  corner  frequency  of  the  Bode  plot  is  the  short 
period  frequency  of  3.0  radians  per  second.  The  response  is 

reduced  due  to  the  effects  of  the  bank  angle.  That  is,  as  the 
aircraft  increases  its  equilibrium  bank  angle,  the  aircraft's 
pitch  rate  response  produces  less  pitch  attitude  change.  This 
can  be  seen  by  examining  the  equation  for  increasing  ^  .  As 

the  aircraft  approaches  90°  of  bank,  pitching  motion  produces  no 
change  in  pitch  attitude. 

The  "’/'Se  is  generally  unchanged  for  differing  values  of  load 
factor  (Figure  B6-B10).  It  contains  a  spike  at  the  same  fre¬ 
quency  as  the  spike  in  the  ©/£<.  transfer  function.  The  short 
period  frequency  is  also  3.0  radians  per  second. 
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The  c,  /  f  response  is  not  altered  s  iqn i  f  i  ran t 1 y  by  load 
factor.  The  response  using  the  uncoupled  lateral  equations  and 
the  coupled  equations  with  zero  bank  angle  is  the  same  (Figure 
B11-B12) .  At  load  factors  greater  than  one,  a  spike  at  the 
phugoid  frequency  occurs  and  the  very  low  frequency  response  does 
not  go  to  zero  as  rapidly  (Figure  B12-B14).  These  two  results 
occur  because  of  coupling  between  axes. 

CONCLUSIONS 

During  steady,  level  turns,  the  frequency  response  between 
NT-33A  aircraft  dynamics  and  control  inputs  is  not  significantly 
altered  when  the  load  factor  is  varied.  In  the  range  of  normal 
input  frequencies,  these  variables  are  especially  consistent. 

Thus,  designs  optimized  for  a  flight  condition  such  as 
straight  and  level  will  retain  the  same  response  characteristics 
in  a  steady,  level  turn,  which  is  a  coupled  flight  condition. 
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APPENDIX  b 

COMPUTER  LISTINGS  AND 
bode  PLOTS 
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!  10 

RcM 

THESIS  AIDt  El  Nib  THE'  'A*  AND  *B*  MATRICES  FOR  THE 

20 

RiM 

BASIC  COUPLED  AIRCRAFT  EQUATIONS. 

i  30 

RiM 

i  40 

HiM 

BY  JIM  PAYNE 

I  50 

RiM 

mahgh/aphil  1982 

60 

RiM 

|  70 

DIM  A(8,8) ,E(2,2) , F(2) , H(2) ,q(2) ,U(8) ,2(2,2) 

80 

DEE 

FN  R(l)  =  INT  (1  *  1000000  +  .5)  /  1000000 

90 

RiM 

100 

HiM 

ENTER  CONSTANTS 

110 

RiM 

120 

HiM 

NT-33A  MACH  =  .4  H  =  SEA  LEVEL 

;  130 

RiM 

140  OK  -  57.29578 
150  G  =  32.1725 
160  B  =  37.54 
170  G  =  6.72 
180  S  =  234.8 
190  WT  =  13700 
200  M  -  WT  /  G 
210  IX  =  238OO 
220  IY  -  21100 
230  12  =  43800 
240  2X  =  480 
250  RiM 

260  HiM  ENTER  VARIABLE  QUANTITIES 
2?0  HMD  HO,  V, DO,  m, AO,  CL, CD 
280  DATA  .002378,447, .0185, 5200,. 9,. 24,. 02 

290  AO  «=  AO  /  DR 
300  «£« 

310  Hi®  EN TER  NON-DIMENSIONAL  DATA 
320  RbH 

330  HEAD  DA,YP,YR,LA,LQ,LE,LV,DV 

340  DATA  .133.0.0.5.9,0,  .>5.  .000183,0 

350  GW  =  CL 

36O  KX  =*  (CD  -  DO)  /  (GL  *  GL) 

370  D3  =>  DA 

380  HEAD  LB.LF.LR.NB.Nt-.NH 

390  DATA  -.09, -.58,. 095, -057. -.008, -.148 

400  HEAD  MD,MQ,NE,MA,MV 

410  DATA  -4. 3, -10. 33, -.90, -.51,0 

420  q(1)  “  Mqiq(2)  =  LU 

430  2(1,1)  =  MA >2(1,2)  =  ME 

440  2(2,1)  =  LAi2(2,2)  =  LE 

450  Hah 

460  RcM  CONVERT  INERTIA  TO  STABILITY  AXIS 
470  HEM 

480  GA  =  GUS  (  -  AO) 

490  SA  -  SIM  (  -  AO) 

500  Il  =  IX*CA*CA+2*2X*OA*SA  +  IZ*SA*SA 

) 
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510  13  =  IZ  *  CA  *  CA  -  2  *  ZX  *  GA  *  SA  +  IX  *  SA  *  SA 
520  ZX  =  (IZ  -  IX)  *  GA  *  SA  +  ZX  *  (CA  *  a  -  SA  *  SA) 

530  IX  -  III 12  =  13 
540  nm 

550  KE-i  CUN  VERT  TO  DIMENSIONAL  DATA 
560  TV  -  0 

570  qq  =  RO  *  V  *  V  *  S  /  2 

58O  DV  =  qq  *  DV 

590  DA  =  qq  *  DA 

600  LV  =  qq  *  LV 

610  LB  -  qq  *  B  *  LB 

620  LA  -  qq  *  LA 

63O  LI  -  qq  *  B  *  LP 

640  Lq  =  qq  *  Iq 

650  LH  =  qq  *  B  *  LR 

660  LE  =  qq  *  LE 

6?0  MV  =  qq  *  C  *  MV 

680  MA  =  qq  *  C  *  MA 

690  MD  =  qq  *  C  *  MD 

700  Mq  =  qq  *  C  *  Mq 

710  ME  -  qq  *  C  *  ME 

720  NB  =>  qq  *  B  *  NB 

730  NB  =  qq  *  B  *  NB 

740  NB  =  qq  *  B  *  NR 

750  YB  =  qq  *  YB 

?60  YR  =  qq  *  YH 

7?0  RiiM 

780  RE-1  ENTER  DIMENSIONAL  DATA 
790  REM 

800  READ  XU,YB,ZU,XE,XW,zW,ZE 

810  DATA  -.0104, -.181  ,-.128,  .620,. 0562,-1. 73,-^. ^ 

820  deft 

83O  HiM  TRANSFORM  BODY  TO  STABILITY  AXIS 
840  REM 

85O  DV  =  -  (XU  *  CA  *  GA  -  (XW  +  ZU)  *  CA  *  SA  +  ZW  *  SA  *  SA 

860  D3  =  -  (Aw  *  CA  *  GA  +  (XU  -  ZW)  *  CA  *  SA  -  ZU  *  SA  *  SA 

870  LV  “  -  (ZU  *  CA  *  CA  -  (ZW  -  XU)  *  CA  *  SA  -  XW  *  SA  *  SA 

880  L3  =  -  (Zw  *  CA  *  CA  +  (ZU  +  XW)  *  CA  *  SA  +  XU  *  SA  *  SA 

890  DE  «  -  (XE  *  CA  -  ZE  *  SA) 

900  L9  =  -  (ZE  *  CA  +  XE  *  SA) 

910  REl 

920  daft  REDEFINE  NASA  DATA  FOR  EQUATIONS  OF  MOTION 
930  REM 

940  YB  -  M  *  V  *  YB 
950  DV  **  DV  *  M 
960  D3  =  D3  *  M  *  V 
970  LV  -  LV  *  M 
980  L3  =  L3  *  M  *  V 
990  DE  **  DE  *  M 
1000  L9  -  L9  *  M 


) 
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1010  MD  -  MI)  *  C  /  (2  *  V) 

1020  hv*  -  Mv*  *  C  /  (2  *  V) 

1030  LD  =  LD  *  C  /  (2  *  V) 

1040  LP  =  LP  *  B  /  (2  *  V) 

1030  LR=LR*B/(2*vj 
1060  NP  =  NP  *  B  /  (2  *  V) 

1070  NR  =  NR  *  B  /  (2  *  V) 

1080  Hid* 

1090  HiM  LOAD  MATRIX  FOR  DELTA  ALPHA  SOVLER 
1100  RiM 

1110  FOR  I  =  1  TO  2 
1120  FOR  J  **  1  TO  2 
1130  E(I,J)  -  Z(I,J) 

1140  NEXT 
115)  NEXT 
1160  RuM 

1170  Ritfl  FIGURE  TRIANGULAR  FACTORS  OF  MATRIX  SOLVER 

1160  am 

1190  NN  -  2 

1200  F(NN)  »  1 

1210  FOR  K  =  1  TO  NN  -  1 

1220  MM  -  K 

1230  FOR  I  =  K  +  1  TO  NN 

1240  IF  ABS  (E(I,K))  <r  =  ABS  (E(MM,K))  THEN  1260 

1230  MM  =  I 

1260  NEXT  I 

1270  F^K)  =■  MM 

1230  IF  m  =*  K  THEN  1300 

1290  F(NN)  =  -  F(NN) 

1300  P  =  e(mm,k) 

1310  E(MM,X)  -  E(K,X) 

1320  E(K,K)  «  P 
1330  FOR  I  =  K  +  1  TO  NN 
1340  E(I,K)  =  -  E(I,K)  /  P 

1330  NEXT  I 

1360  FOR  J  =  K  +  1  TO  NN 
1370  T  =  E(MM,J) 

1380  E(MM, j)  -  E(K,J) 

1390  E(K,J)  -  T 

1400  IF  T  =  0  GOTO  1440 

1410  FOR  I  =  K  +  1  TO  NN 

1420  E(I,J)  -  E(J.,J)  +  E(I,K)  *  T 

1430  NEXT  I 

1440  NEXT  J 

1430  NEXT  K 

1460  R*M 

1470  HFK  THE  MAIN  LOOP 
1480  RfcM 
1490  READ  N 
1500  DATA  4 

) 
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1510  FK  =  ATN  (  S*R  (N  *  N  -  1))  +  SGN  (N)  -  1)  *  1 . 5707963 

1520  M  =  G  *  TAN  U'*0  /  V 

1530  FF  =  *  *  C  *  SIN  (FE)  /  (2  *  V) 

1540  H(l)  =  -  *(1J  *  *’*’ 

1330  H(2)  =  -  *(2)  *  KF  +  (N  -  1)  *  GW 

I36O  14E  =  *1  *  SIh  (FE) 

13?0  RE  -  W  *  COS  (FK) 

1390  RaM 

1590  tm  back,  substitute  to  find  delta  alpha 

1600  Kl<« 

1610  FOR  K  **  1  TO  NN  -  1 
1620  MM  =  F(K) 

1630  SB  =  H(MM ) 

1640  H(MM)  =  H(K) 

1650  H(k)  =  Si 

1660  FOR  I  =  K  +  1  TO  NN 

1670  H(l)  =  H(l)  +  E(l,K)  *  SS 

1680  NEXT 

1690  NEXT 

1700  FOR  J  =  1  TO  NN  -  1 
1710  K  ■»  NN  -  J  +  1 
1720  H(K)  =  H(K)  /  E(K,K) 

1730  ss  =  -  h(k) 

1740  FOR  I  =  1  TO  NN  -  J 
1730  H(l)  -  H(l)  +  E(I,K)  *  SS 
1760  NEXT 
1770  NKXT 

1780  H(l)  =  H(l)  /  E(l,l) 

1790  AD  -  H(,l) 

1800  REN  CHECK  THAT  DRAG  THRUST 
1810  L  =■  N  *  WT 

1820  CL  -  2  *  L  /  (RO  *  V  *  V  *  S) 

1830  D=RO*V*V*S*  (D0+XK*CL*CL)/2 

1840  IF  D  >  TM  THEN  PRINT  1  PRI  "N  =  "jNs",  DRAG  >  THRUST"»  END 

I83O  RhM 

1860  REM  MOVE  STABILITY  AXIS  BY  DELTA  ALPHA 
I87O  Re*1 

1880  CA  -  COS  (  -  AD) 

1890  SA  -  SIN  (  -  AD) 

1900  Dl  -  DV  *  CA  *  CA  -  (LV  +  DA  /  V)  *  SA  *  CA  +  (LA  /  V)  *  SA  *  SA 

1910  D3  -  (DA  /  V)  *  CA  *  CA  +  (DV  -  LA  /  V)  *  SA  *  CA  -  LV  *  SA  *  SA 

1920  D3  -  D3  *  V 

1930  Y4  -  Yp  *  CA  -  YR  *  SA 

1940  Y6  -  YR  *  CA  +  YP  *  Sa 

1930  LI  -  LV  *  CA  *  CA  -  ((LA  /  V)  -  DV)  *  SA  *  CA  -  (DA  /  V)  *  SA  *  SA 

1960  LX  =  (LD  /  V)  *  CA  *  CA  +  (DD  /  V)  *  SA  *  CA 

1970  LX  -  LX  *  V 

1980  L5-W*CA+D«*SA 

1990  L3  -  (LA  /  V)  *  CA  *  CA  +  (LV  +  DA  /  V)  *  SA  *  CA  +  DV  *  SA  *  SA 
2000  L3  -  L3  *  V 


) 
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2010  14  =  (IZ  -  IX)  *  CA  *  SA  +■  ZX  *  (CA  *  a  -  SA  *  SA) 

2020  11  =  lA*CA*CA+2*ZX*SA*CA+IZ*SA*SA 
2030  L2  =  (LB  /  V)  *  CA  -  (NB  /  V)  *  SA 
2040  L2  =  L2  *  V 

20^0  L4  =  Lp*CA*CA-  (LR  +  NP)  *  SA  *  CA  +  NR  «  5A  *  SA 
2060  I3  =  IZ*CA*CA-2*ZX*SA*CA  +  IX*SA*SA 
2070  L6«LR*CA*CA-  (NR  -  LP)  *  SA  *  CA  -  NP  *  SA  *  SA 
2080  N2  =  (NB  /  V)  *  CA  +  (LB  /  V)  *  SA 
2090  N2  =  N2  *  V 

2100  N4  =  NP*CA*CA-  (NR  -  LP)  *SA*CA-LR*SA*SA 

2110  N6  =  NR  *  CA  *  CA  +  (LR  +  NP)  *  SA  *  CA  +  LP  *  SA  *  SA 

2120  Ml  -  MV  *  CA  -  (MA  /  V)  *  SA 

2130  M4  =  (MD  /  V)  *  CA 

2140  M4  -  M4  *  V 

2150  M3  -  (MA  /  V)  *  CA  +  MV  *  SA 

2160  M3  -  M3  *  V 

2170  D9=’DE*CA-LE*SA 
2180  L9  =  LE  *  CA  +  DR  *  SA 
2190  RbJM 

2200  HUM  FIGURE  ELfcMENTS  OF  'A*  MATRIX 
2210  HiiM 

2220  A(l,l)  -  (TV  *  COS  (AO  +  AD)  -  Dl)  /  M 
2230  A(l,2)  =  RE  *  V 

2240  All, 3)  =  (  -  D  *  SIN  (AO  +  AD)  +  L  -  M  *  qE  *  V  -  D3)  /  M 

2250  A(l,8)  =  -  G 

2260  A(2,l)  =  HE  /  V 

2270  a(2,2)  =  YB  /  (M  *  V) 

2280  A  (2,4)  =  Y4  /  (M  *  V) 

2290  A (2, 6)  -  -  1  +  Y6  /  (K  *  V) 

2300  A (2, 7)  “  G  *  COS  (FE)  /  V 

2310  A(3,l)  -  (  -  TV  *  SIN  (AO  +  AD)  -  LI  -  (4E  *  M)  /  (M  *  V  +  LX) 

2320  A (3, 3)  -  (  -  L3  -  D  *  COS  (AO  +  AD)  -  D)  /  (M  *  V  +  LX) 

2330  A (3, 5)  =  (V  *  M  -  L5)  /  (M  *  V  +  LX) 

2340  A (3, 7)  -  (  -  WT  *  SIN  (FE))  /  (M  *  V  +  LX) 

2350  A (4, 2)  =  L2  /  II 

2360  A (4,4)  -  (IA  +  14  *  qE)  /  II 

23?0  A (4, 5)  =  ((IY  -  13)  *  RE)  /  II 

2380  A(4,6)  -  (L6  +  (IY  -  13)  *  qE)  /  II 
2390  I  -  14  *  (14  /  II) 

2400  a(6,2)  -  (N2  +  14  *  A(4,2))  /  (13  -  I) 

2410  A(6,4)  -  (N4  +  (II  -  IY)  *  qE  +  14  *  A (4,4))  /  (13  -  l) 

2420  A(6,5)  =  (  -  14  *  RE  +  14  *  A(4,5))  /  (13  -  i) 

2430  a(6,6)  **  (N6  -  14  *  qE  +  14  *  a(4,6))  /  (13  -  l) 

2440  I  =  1  /  II 

2430  A (4, 2)  =  I  *  A(6,2)  +  A(4 ,2) 

2460  A(4,4)  -  I  *  A(6,4)  -t-  A(4,4) 

2470  A (4, 5)  -  I  *  A  (6, 5)  +  A  (4, 5) 

2480  A(4,6)  =  I  *  A(6,6)  +  A(4,6) 

2490  A(5,l)  =■  (Ml  +  M4  *  A(3,l))  /  IY 

2300  A  (3, 3)  =  (M3  +  M4  *  A  (3, 3))  /  IY 
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2510  A(,5.*0  =  Ul3  -  II)  *  Hi,')  /  IY 
2520  A 15, 5)  =  (M*  +  M4  *  A (3, 5))  /  IY 
2530  A15.6)  =  2  *  14  *  RE  /  IY 
2540  A(5,7)  =»  M4  *  A (3, 7)  /  IY 
255)  A(7,4)  =  1 

2560  A(7,8)  =  OK  *  SIN  (FB)  +  RE  *  CCS  (FE) 

2570  A (8, 5)  =  COS  .FE) 

2580  A(8,6j  =  -  SIN  (FB) 

2590  A(8,7)  =  -  QE  *  SIN  (FE)  -  HE  *  COS  (FE) 

2600  REM 

2610  Hi*  FIGURE  ELEMENTS  OF  'B1  MATRIX 

2620  REM 

2630  U(l)  =  -  D9  /  M 

2640  U(3)  =  -  L9  /  (M  *  V  +  LX) 

2650  U(5)  =  (ME  +  M4  *  U(3))  /  IY 
2660  REM 

2670  REM  PRINT  RESULTS 

2680  REM 

2690  PRINT  "THE  'A*  MATRIX  EY)R  LOAD  FACTOR,  N  =  "iN 
2700  PRINT  «  PRINT  "BANK  ANGLE,  FE  =  “;FE  *  DR; "  DEGREES" 

2710  PRINT  "DELTA  ALPHA,  AD  =  ";AD  *  DR;"  DEGREES" 

2720  PRINT 

2730  FOR  I  =  1  TO  8 

2740  FOR  J  =  1  TO  8 

2750  PRINT  FN  H(A(l, J)); "  "; 

2760  NEXT 
2770  PRINT 
2780  NEXT 
2790  PRINT 

2800  PRINT  "THE  'B*  MATRIX  TRANSPOSED  IS" 

2810  PRINT 

2820  FOR  I  =  1  TO  8 

2830  PRINT  FN  R(U(l)) ; "  "; 

2840  NEXT 

2850  PRINT  :  PRINT 

) 
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HUN 

THE  'A*  MATRIX  FOR  LOAD  FACTOR,  N  =  4 

BANK  ANGUS,  FE  =  75.5224885  DEGRESS 
DELTA  ALPHA,  AD  =>  7.31234-296  DEGREES 

-.064253  31.150 889  -IO5.68743I  0000  -32.I725 
1.56E-04  -.181  000-1  .017994  0 

-1.428E-03  0  -1.702174  010  -.069689  0 

0  -7.134195  0  -2.07948  -.064927  .289606  0  0 

-1.766E-03  0  -8.017326  .064257  -1.94903  -.OI5636  .039692 

0  3.678639  0  .22388  7.36E-03  -.338186  0  0 
0001000  .278755 
0  0  0  0  .25  -.968246  -.2?8755  0 

THE  'B'  MATRIX  TRANSPOSED  IS 

-5.829034  0  -.10026  0  -15.93 2007  000 
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FREQUENCY  (RflO/SEC) 


THETfl  /  DELTA  E  FOR  NT-33R 


FREQUENCY  (RflO/SEC) 


THETA  /  DELTA  E  FOR  NT-33A, 


FREQUENCY  C  RRO/SEC ) 


THETR  /  OELTft  E  FOR  NT-33R, 


FREQUENCY  (RfiO/SEC) 


THETfl  /  OELTfl  E  FOR  NT-33A, 


FREQUENCY  ( RflO/SEC ) 


FREQUENCY  tRRD/SEC) 


FREQUENCY  IRflO/SEC) 


RflO/SEC ) 


FREQUENCY  (RflO/SEC) 


Figure  Bll. 


FREQUENCY  IRflO/SEC) 


FREQUENCY  IRflO/SEC) 


APPENDIX  C 
DATA  PLOTS 
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COOPER  - HARPER  VS  OMEGA  A!'  1G 
SRFTD  TV  TRACKING  TASK 
N  /  ALPHA  =  29  G  /  RAD 
0.7  DAMPING  RATIO  .  1G/2ND  ORDER 

(  PILOT  ffl  RATIN6S  )  - 


ratings  c-h  ratings 


RATINGS  C-H  RATINGS 


COOPER-HARPER  VS  GMEGR  RT  2G 
SRFTD  TV  TRACKING  TASK 
N  /  ALPHA  =  29  G  /  RAD 
0.7  DAMPING  RATIO,  2G/2ND  ORDER 


OMEGA  NATURRL  FREQUENCY  (rad/sec) 

Figure  ce .  c  0  0  P  ER  -  HR  R  P  E  R  VS  OMEGA 


74 


RATINGS  C-H  RATINGS 


COOPER-HARPER  VS  OMEGA  RT  2G 
SRFTD  TV  TRACKING  TASK 
N  /  ALPHA  =  29  G  /  RAD 
0.7  DAMPING  RATIO,  2G/2ND  ORDER 

(  PILOT  #3  RATINGS  ) 


NOTE: 

scrips  indicate  the 
r  of  overplotted 

Super 

numbe 

<1 

_ 

Pjjl| 

|| 

i 

i 

2 

*  J 

t 

1 

l 

i 

c_> 


L 


1  1 

(  2nd  ORDER  CURVE  FIT  ) 

r 

< 

: 

V 

• 

) 

:  : 

0  2  4  6  3  18  12 


OMEGA  (fcy,  NATURAL  FREQUENCY  (rad/sec) 
Figure  ci.  COOPER -HRRPER  VS  OMEGA 
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RATINGS  C-H  RATINGS 


COOPER-HARPER  VS  OMEGA  RT  2G 
SflFTD  TV  TRACKING  TRSK 
N  /  ALPHA  =  29  G  /  RAD 
0.7  DAMPING  RATIO,  2G/2ND  ORDER 


(  COMPOSITE  ) 


LEGEND 

1 . 

0:  prior 

1 

2. 

0:  PILOT 

£ 

3. 

A:  PILOT 

4. 

* :  PILOT 

1 

s . 

Z  :  PILOT 

~) 

6 . 

X;  PILOT 

0  2  4  6  8  10  12 


OMEGA  <42  >  NATURAL  FREQUENCY  (^d/sec) 
Figure  cs.  C 0 0 P E R - H A RPE R  VS  OMEGA 
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-H  HOT INGS  C-H  RATINGS 


Ratings  c-h  ratings 


RATINGS  C-H  RATINGS 


COOPtR-MARPER  VS 
NT-33R  PROJECT 
N  /  ALPHA  = 
0.7  DAMPING  RATIO 


S  OMEGA  AT  2G 
HUD  TRACKING 
29  G  /  RAD 
,  2G/4TH  ORDER 


PILOT  #3  RATINGS 


Superscrips  indicate  the 
number  of  overplotted 
data  pionts. 


RMS  ERROR  (X  1000) 


«««»  -  Jr,"5, 


F/e  w/a 


RMS  ERROR  (X  1000) 


PILOT  *5  / 
NT-53P  PROJECT  I 
N/RLPMR  : 
0.7  ORMP 


2  3  4 

COOPER-HRR 


urncn 


RMS  ERROR  <X  1000) 


1 


PILOT  #1  /  SRF7D  1G 
NT-33R  PROJECT  HUD  TRRCKI 
N/RLPHA  =  29  G/RRD 
0.7  DAMPING  RATIO 


NG 


f 


NATURAL  FREQUENCY  (RAD/SEC 


Figure  C20.  RMS  VERSUS  SHORT  PERIOD  NATUR 


RMS  ERROR  (X  1000) 


RMS  ERROR  <X  1000) 


PILOT  #3  /  5RFTD  1G  DRTR 
NT-33R  PROJECT  HUD  TRACKING  TRSK 
N/flLPHR  =  29  G/RRD 
0.7  ORMPING  RATIO 


LEGEND 


10  12  U 


NRTURRL  FREQUENCY  (RRD/SEC ) 


Figure  C22.  RMS  VERSUS  SHORT  PER  I OD~ NR TURRL  TRFOUENC Y 


rsjco 


700 


Figure  C24.  RMS  VERSUS  COOPER  HARPER  RPT  I  NG/NATURflL  FREQUE 


3AFTD  2G  DAT 
NT -3:-5fi  PROJECT  HUD  TRACK 
N/ALPHA  =  29  G/RRD  FS 

0.7  DAMPING  Rfl 
legend 


RATING  NAT  FREQ  (RRD/3EC) 


ERROR  COUNT 


PILOT  #Z  /  NT-33R  1G  DRTR 
NT-33R  PROJECT  HUD  TRRCKING  TRSK 
N/RLPHR  =  29  G/RRD 
0.7  DAMPING  RATIO 


300 

275 

250 

225 

200 

175 

150 

125 

100 

1 

f 

Figure  c27. error  COUNT  VERSUS  COOPER-HRRPER  RRTING 
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COOPER-HRRPER  RATING 


ERROR  COUNT 


NT-33R  1G  DRTR 

IUD  TRRCKING  TASK 

:  29  G/RRD 

NG  RRTIO 

1  FGFNP. 

2  3 

COOP 


Figure  C28.  £RR 


ERROR  COUNT 


» 
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PILOT  #2  /  NT-33R  1G  DPTfl 
NT-33R  PROJECT  HUD  TRACKING  TRSK 
N/RLPHR  =  29  G/RRD 
0.7  OflMPING  RRTIO 

300 

275 

250 

22S 

200 

175 

150 

125 

100 

0 

Figur*  C29. 

ERROR  COUNT  VERSUS  SHORT  PER IOO  NRTURRL  FREQ 
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NATURAL  FREQUENCY  (RAO/SEC) 


NT-33R  1G  DATA 
UD  TRACKING  TR5K 
29  G/RRD 
NG  RRTIO 


Table  D1 

SAFTD  STRIP  CHART  OUTPUTS 


Parameter 

Description 

Strip  Chart 

Range 

Fs 

Long.  Stick  Force 

+  30  lbs 

Fs/Nz 

Stick  Force  per  g 

0  to  10  lbs  per  g 

Nz 

Normal  Acceleration 

0  to  4  g  ’  s 

RMS 

Root  Mean  Square 

0  to  2 

Test  Pulse 

SAFTD  Tracking  Task 

event 

Ae 

Error  Between 

+  10  deg 

Tracking  Task  and 

Pitch  Angle 

e 

Pitch  Angle 

+  10  deg 

oC 

Angle  of  Attack 

+  10  deg 

Note:  Strip  charts  were  played  back  at  0.5mm  per  second  at  a 


scale  factor  of  0.2  volts  per  division. 


Table  D2 

NT-33A  STRIP  CHART  OUTPUTS 


VSS  Scale 


Digital 

Channel  No. 

Parameter 

Description 

Output 

P.P. 

Factor 
(x  per  volt) 

3* 

Vi 

Indicated  Velocity 

106 

27.5  Kt 

18 

Fes 

Long.  Stick  Force 

26 

20.0  lb 

2 

Nz 

Normal  Acceleration 

16 

0.5  g 

24 

«v 

AOA  (vane) 

38 

2 . 0  deg 

4 

q 

Pitch  Rate 

4 

5.0  deg/sec 

5 

sin© 

Pitch  Angle 

28 

5.88  deg 

11 

£ec 

Tracking  Task  Error 

D/A#l 

2.5  deg 

1** 

hp 

Pressure  Altitude 

107 

1 , 300  feet 

Notes:  1.  *  Zero  knots  equalled  10.00  volts. 

2.  **  720  feet  equalled  8.45  volts. 

3.  Digital  channels  one  through  12  had  adjustable  gain  xl,  x2, 
x5  and  xlO.  Gains  were  xl . 

4.  Aircraft  data  tapes  were  played  back  at  1.0mm  per  second. 
Strip  chart  sensitivity  was  0.1  volt  per  division. 
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Table  El 


ft 

PILOT  BACKGROUND 


PAYNE,  JAMES  M.,  Captain,  USAF 

EDUCATION:  BS  -  General  Engineering,  US  Air  Force  Academy,  CO 


FLIGHT  TIME:  F-4C/E 

602 

F-5E 

521 

Other 

870 

Total 

1993 

OPERATIONAL  EXPERIENCE:  F-4E  pilot  for  one  year  in  the  334th 
Tactical  Fighter  Squadron  at  Seymour  Johnson  AFB,  NC;  F-4C/E 
pilot  for  fifteen  months  in  the  57th  Fighter  Interceptor 
Squadron,  Keflavik  NS,  Iceland;  F-5E  Aggressor  Dissimilar  Air 
Combat  Tactics  Instructor  Pilot  for  two  years  in  the  527th 
Tactical  Fighter  Training  Aggressor  Squadron,  RAF  Alconbury, 
United  Kingdom. 


'  ..  3 
•  .  1 


STEWART,  CHARLES  R. ,  Captain,  USAF 

EDUCATION:  BS/MS  -  Aeronautical  and  Astronautical  Engineering 

Purdue  University,  West  Lafayette,  Indiana 

FLIGHT  TIME:  A-10  1325 

Other  331 

Total  1656 

OPERATIONAL  EXPERIENCE:  A-10  pilot  for  two  years  in  the  356th 
Tactical  Fighter  Squadron  at  Myrtle  Beach  AFB,  SC  and  then  served 
another  two  year  tour  as  an  A-10  Instructor  Pilot  with  the  92nd 
Tactical  Fighter  Squadron  at  RAF  Bentwaters,  United  Kingdom. 


TOMENY,  TERRY  E. ,  Captain,  USAF 

EDUCATION:  BS/ME  -  Mechanical  Engineering, 

Rensselaer  Polytechnic  Institute,  Troy,  NY 

FLIGHT  TIME:  T-38A  1350 

RF-4C  1000 

Other _ 100 

Total  2450 

OPERATIONAL  EXPERIENCE:  T-38A  Instructor  Pilot  for  four  years  at 
Moody,  AFB,  GA.  and  Vance  AFB,  Okla;  RF-4C  Instructor 
Pilot/Flight  Examiner  for  four  years  at  Zweibrucken  AB,  Germany. 
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VITA 


James  M.  Payne  was  born  on  July  2nd,  1952,  in  Eugene, 

Oregon,  and  grew  up  in  Hastings,  Michigan.  Graduation  from 
Hastings  High  School  was  followed  by  undergraduate  schooling  at 
the  United  States  Air  Force  Academy.  He  graduated  from  the 
Academy  in  1974  with  a  Bachelor  of  Science  degree  in  General 
Engineering  and  a  commission  as  a  Second  Lieutenant  in  the  United 
States  Air  Force.  He  completed  Undergraduate  Pilot  Training  at 
Laughlin  Air  Force  Base,  Texas,  in  1975,  and  F-4C  combat  training 
at  Luke  Air  Force  Base,  Arizona,  in  1976.  After  one  year  of 
flying  F-4E's  at  Seymour  Johnson  Air  Force  Base,  North  Carolina, 
he  was  sent  to  Keflavik  Naval  Station,  Iceland,  for  a  fifteen 
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